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Abstract 

We study the effective action for gravity obtained after the integra- 
tion of scalar matter fields, using the local momentum representation 
based on the Riemann normal coordinates expansion. By considering 
this expansion around different space-time points, we also compute 
the non-local terms together with the more usual divergent ones. We 
discuss the applicability of our results to the calculation of particle pro- 
duction rates in cosmological backgrounds and compare this method 
with the traditional Bogolyubov transformations. 
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1 Introduction 

One of the most interesting aspects of quantum field theory in the presence of 
a external classical background is the possibility of particle creation. When a 
external field J is present, the vacuum state (with zero particle number) can 
become unstable and eventually decay by particle emission. In other words, 
if the system is prepared in the vacuum state in the remote past |0,m), the 
probability that the system remains in the same state in the far future will 
be in general different from one: |(0, m|0, out)j\^ ^ 1. 

The first works about this subject can be traced back to the early fifties 
in the context of quantum electrodynamics (QED) In particular, it was 
shown that intense electric fields could give rise to the creation of electron- 
positron pairs. The corresponding particle creation rates were obtained by 
means of effective action (EA) techniques. Within this formalism the vacuum 
persistence amplitude (0, m|0, owt) j is given by means of the path integral. 
In the simple case of a real scalar field in a external gravitational background, 
the EA is given by: 



where Oxyijn?) = {—Oy — m'^ — ^R{y) + ie)6'^{x,y) with 6^{x,y) the covariant 
Dirac delta. As usual we introduce the non- minimal coupling ^R(f)^ in order 
to include the case with conformal invariance ^ = 1/6. As we see from this 
expression, the functional integration is performed only on the scalar fields, 
i.e, they will appear inside the loops, whereas the gravitational field g^u is 
considered as classical and hence it will only appear in the external lines. 
As a consequence, TVfs'^i/] is the generating functional of Green functions 
with external gravitational lines and loops only of matter fields. Another 
interesting point is that as far as the classical action is quadratic in the 
scalar fields, the one loop calculation is exact, i.e there is no higher loops 
contribution to the EA. 

It is now relatively easy to compute the pair production probability Pea 
that according to the above discussion will be nothing but: Pea = 1 — 
\{0,in\0,out)g^J'^. From eq.|l| assuming W[gf^u] to be small, we obtain: 




(1) 



Pea = 1 - e 



21mW[g,,,] 



(2) 
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Thus the imaginary part of the effective action provides the particle creation 
probabihties. 

Since the first developments of quantum field theory in curved space- 
time, the topic of particle creation by gravitational fields has received much 
attention. The pioneering works were due to Parker who studied particle 
production in Robertson- Walker geometries by means of the mode-mixing 
Bogolyubov method f^. Some other important works are those of Zel'dovich 
and Starobinsky 0, ^ that extended this technique to homogeneous but 
anisotropic cosmological space-times. They show that particle production 
provides an effective mechanism for anisotropy damping in the early universe. 
But probably it is Hawking radiation the best known example of particle 
production from a gravitational field. Most of these calculation were done by 
means of Bogolyubov method. However in this work we present an alternative 
technique for the calculation of particle creation based on the gravitational 
EA. 

The gravitational effective action eq.|l| is in general a non-local functional 
on the metric tensor and the connection. Although the divergent local terms 
have been extensively considered in the literature Q, the finite non-local 
contributions have been studied only more recently [0, ^ 0; il- The non- 
local terms are responsible for the imaginary part and therefore also for the 
particle creation. The calculation of the EA is in general very difficult for an 
arbitrary space-time geometry and accordingly we have to rely on some kind 
of approximation. The most traditional one is the well-known perturbation 
theory in which the metric tensor is splitted as follows: g^y = rj^^ + nh^^. The 
EA is then expanded in powers of the coupling k. This method however is not 
generally covariant. We will present an alternative (covariant) approximate 
method in which the EA is expanded in powers of curvature tensors. The 
curvature expansion will be generated in turn from the Riemann normal 
coordinates expansion. This kind of curvature expansions was first considered 
by Barvinsky and Vilkovisky 

2 Non-local gravitational effective action 

Let us consider a smooth manifold in which geodesies do not intersect so that 
we can perform the well-known Riemann normal coordinates expansion 
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for the metric tensor: 



+ 



I 2 



V/ + Oid') (3) 



where is the origin of normal coordinates. When substituting this expan- 
sion in eq.[^ we will generate the curvature expansion we are looking for. The 
operator Oxy defined above can then be splitted in a free part: 



Axy{m^) = {-al -m^ + te)6%x, y) 



(4) 



with Dg = ri'^'^d^dy, and an interaction term that contains all the curvature 
dependence: 



-^R\{y,)yPdl + ^R^/^{yo)yVdf,dy - ^R{yo) 
1^ 1^ 1 ^ 1 ^ 



20 



20 



20 



-i?„A(yo)i?V7(2/o) + ^Rj'pxiyo)R\%iyo) + i^Ra%xiyo)R\%{yo) 



45 

+ ^R^f,,,''{yo) + ^R.0-,{yo) ) y'^yS'dl 
1 
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R^p\-M^yo) 



+ ^R;Ay^)R\\{yo)) y'y'y'y^'didi - ^-R.Myo)y''y^ 

Accordingly we can write: 

= ^Tr log 0(m2) = ^Trlog(A + B) 



5%x,y) + ... 

(5) 



(6) 



It is then possible to expand in a power series in the interaction operator 
B. With that purpose we first take the derivative of the EA with respect 
to the mass parameter, in order to avoid the problem of series expansions of 
logarithms of operators [§]: 



d 



dm? 



rp 

2 ^0(m2) 



-^Tt{{A + B)-'] 



— Tt{A-^ - A-^BA-^ + A-^BA-^BA-^ - ...) (7) 
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As far as the B operator is at least linear in the curvatures and we are only 
interested in the expansion up to order quadratic in the curvatures, we will 
only have to keep the first three terms. We now face the main problem of 
normal coordinates expansion. As we can see from eq.^ all the curvature 
tensor are evaluated at the same point i/q, and therefore eq.^ can hardly 
gives rise to a non-local expression. However, it is possible to circumvent 
this problem by means of a point splitting procedure in the quadratic terms. 
Using again the Riemann normal coordinate expansion we can write: 

n{yo)n{yo) = n{yo)n{y) + (niyo)Vn{yo)y + n{yo)V^n{yo)yy + ...) 

+ (n{yo)'JZ{yo)'JZ{yo)yy + 'JZ{yo)V'JZ{yo)'JZ{yo)yyy +...) + ... (8) 

where TZ collectively denotes the scalar curvature or the Riemann and Ricci 
tensors. There are several ways of performing the splitting (according to 
the chosen y point), however they all will differ in higher order terms in 
derivatives. 

Evaluating the functional traces in eq.|^ by using dimensional regulariza- 
tion, we obtain the final result for the gravitational EA up to 0{TZ'^) for a 
real scalar field. We show the masless limit [^: 

+ R{x)T{n)R{x)^ + o{n') (9) 

with r(n) = A^, - log(n//i2) and A", = 2/e + log47r - 7 the usual way of 
parametrizing the poles in dimensional regularization. The action of the r(n) 
operators should be understood through the corresponding Fourier transform 
in normal coordinates, i.e: 

log (H) nv.) = /<''-(g,^->og (^) n^) (10) 

We see that the massless limit is regular ^ and in this case, the above 
expression for the EA is exact up to quadratic order. 



5 



3 Particle production in cosmological back- 
grounds 

In this section we are interested in calculating the imaginary part of the 
gravitational EA in eq.^that, as shown in the introduction, provides us with 
the particle production probability. 

Let us consider a general Bianchi I type metric: 

ds^ = a^{T)dT'^ - a^{T)gij{T)dx'dx^ (11) 

As far as this metric is homogeneous, it is possible to write the different 
curvature tensor in such a way that they only depend on the time coordinate, 
and therefore: 

Im / dS^^e^^-ny,) log (^^^) ^(^) (12) 

= Im I rfx°^e^P«-"7^(l/o) log (^^^) ^(^°) = -^7^(l/o)7^(l/o) 

Using this result it is possible to write the imaginary part as a purely local 
expression: 

Here C^^pcr is the Weyl tensor and we have discarded a Gauss-Bonnet term. 
In order to do so, it is necessary that space-time is asymptotically flat or at 
least that asymptotically (r — *• ±cxd), the scale factor behaves as: a(r) ~ \t\" 
with a > —1/2. The integrand can be readily interpreted as the particle pro- 
duction probability per unit time and volume. Therefore, unlike Bogolyubov 
method that only provides the spectrum of the particle produced, the EA 
method allows us to know the instantaneous particle production rate. This 
rate will give rise to the same total amount of particle as the Bogolyubov 
method in contrast with other methods such as Hamiltonian diagonalization 



(see [§] and references therein) which produce vastly more particles. Eq.l3 
is just the general form of the Zel'dovich and Starobinski |^ result. When 
the metric tensor reduces to the Friedmann-Robertson- Walker (FRW) form 
(with vanishing Weyl tensor), the probability will only depend on the scalar 
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curvature. It is then easy to see that in the case of conformal invariance 
= 1/6) there is no particle production, as well as in a radiation dominated 
universe for which i? = 0. 

It is possible to compare (for specific FRW models with ^ = and complex 
scalar field), the particle production probabilities per unit volume {pea and 
Pbog) derived with both methods: 

/d'^ k 1 /" 

^^log(l+ < Nk >) a.nd Pea ^ ^ J dTa\T)R\T) (14) 



As an example let us take the model considered in |jTO[ with a^(r) = A^t"^ + 
, where A and B are arbitrary constants. The number density of cre- 
ated pairs in the k mode is given by the Bogolyubov method < Nj. >= 
exp(— '^^ ^ ). From eq.l4 we extract the following analytical results: 



PBGC = PEA = (15) 

i.e. both methods give the same probabilities. Comparison with other mod- 
els can be found in [|ll| . Finally it is worth mentioning that the EA method 



also provides the particles spectra |TTI] and could be useful in those con 



texts in which Bogolyubov transformations have been traditionally applied; 
such as graviton production or the generation of anisotropics in inflationary 
cosmological models |T2| . 

Acknowledgements: This work has been partially supported by Min- 
isterio de Educacion y Ciencia (Spain) CICYT (AEN96-1634). 



7 



References 

[1] J. Schwinger, Phys. Rev 82, 664 (1951) 

[2] L. Parker, Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. 183, 1057 (1969) 

[3] N.D. Birrell and P.C.W. Davies Quantum Fields in Curved Space, Cam- 
bridge University Press (1982) 

[4] Y.B. Zel'dovich and A. A. Starobinski, Pis'ma Zh. Eksp. Tear. Fiz 26, 
373 (1977) {JETP Lett. 26, 252 (1977)) 

[5] J.B. Hartle and B.L. Hu, Phys. Rev. D20, 1772 (1979) 

[6] A.O. Barvinsky and G.A. Vilkovisky, Nucl. Phys. B282 (1987), 163 
(1987); B333, 512 (1990) 

A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky and V.V. Zhytnikov, 
Nucl. Phys. B439 (1995) 

[7] I.G. Avramidi, Nucl. Phys. B355, 712 (1991) 

[8] A. Dobado and A.L. Maroto, preprint |hep-th/ 97121981 

[9] A.Z. Petrov, Einstein Spaces, Pergamon, Oxford (1969) 

[10] G. Schafer, J. Phys. A12, 2437 (1979) 

[11] A. Dobado and A.L. Maroto, preprint |gr-qc/9803076 

[12] V.F. Mukhanov, H.A. Feldman and R.H. Brandenberger, Phys. Rep. 
215, 203 (1992); Y. Shtanov, J. Traschen and R. Branderberger, Phys. 
Rev. D51, 5438 (1995) 



8 



